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ABSTRACT 



The classical working states of the airscrew were 
studied employing the three-dimensional flow visualization 
tunnel at the Naval Postgraduate School, Monterey, Cali- 
fornia. The various states were simulated in an attempt to 
gain insight into some of the present problems in heli- 
copter rotor design. This method was extended to include 
two particular phenomena, the critical part of the vortex 
ring state and Mangier' s distribution of induced velocity 
about a rotor in the helicopter mode. Using a feasible 
distribution ascertained from the flow visualization, the 
inducea velocity disi r. ibuLiun was calculated fer several 
points in the vortex ring state. 
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TABLE OF SYMBOLS 



Specified values are used in subsequent calculations. 

steady velocity of ascent or descent 
V induced velocity at blade element dr 

r radical distance of element dr from rotor axis 

R rotor radius = 30 feet 

B number of blades 

C chord of blade (considered constant) 

0 geometric blade angle measured from no-lift line 

= .07 radians 

a angle of attack measured from no-lift line 

H. a - 8 or u/o:r 

0) angular velocity of blades, considered constant = 

21 radians/sec 

cr solidity or ratio of blade area to disk area 

BC/ttR = .05 

p density of air 

elemental coefficient of lift 

T thrust along the shaft 

dT elemental thrust due to B blade elements at radial 

distance r 

X non-dimensional rotor radius r/R 

A slope of lift curve = 3 

u flow through rotor; ” v for upwash; v - for 

downwash 
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I. 



INTRODUCTION 



The last decade has witnessed the emergence of the 
helicopter's increasingly important role in military and 
civil aviation. Militarily, it has become the backbone 
of the Army's tactical doctrine. The other services have 
also found it increasingly reliable in their own particular 
operations. Civil aviation has barely scratched the surface 
of the potential applications of the helicopter. It has 
been used as a shuttle service between large urban areas 
and their major airports. Recent studies have shown that 
this service could be profitably extended to short and 
medium range eperations, in compet tv’ii-h rhf^ paKRenger 
jets. Like the proverbial tortoise and the hare, the heli- 
copter, which can transport the passenger from door-step 
to door-step, is often a faster mode of transportation than 
the jet liner. 

The demand for faster and more versatile helicopters 
has renewed interest in low speed aerodynamics after a long 
emphasis on supersonic, hypersonic, and space applications. 

The flow field of a helicopter rotor is extremely diffi- 
cult to model mathematically. The periodicity of the rotor, 
areas of reverse flow, and the non-uniformity of the inflow 
are only a few of the problems which are confronted when 
tryirig to calculate flow and performance parameters. 
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One tool used to aid in these calculations is flow 



visualization. From a visual study of the flow in a rotor 
system, a feasible mathematical expression for the different 
variables may be determined. In this paper flow visuali- 
zation is used to examine the working states of the rotor 
and to assist in the calculation of induced velocity in 
the vortex ring state, that is to say, recirculatory flow 
at negative speeds of advance. 
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II. FLOW VISUALIZATION 



With an increasing interest in STOL/VSTOL aircraft, 
some of the unsolved low speed aerodynamic problems, which 
have fallen by the wayside in the space age, are receiving 
more attention. As discussed previously, one of the more 
useful tools in selecting a mathematical model is flow 
visualization. With this in mind a flow visualization 
facility has recently been completed at the Naval Post- 
graduate School. This facility is modeled after the one 
described in Reference 5. 

Several methods were attempted to obtain the cloud used 
as the visuaji 7:atir>Ti inr^ri-inm. effective V7as found 

to be the vaporization of silicone-based turbine solvent. 

This fluid produces a dense white coherent cloud. Low- 
pressure air forces the fluid to drop on a heated surface 
in a locally manufactured cloud generator heated by a 
2400-watt hot plate. As the vaporized fluid collects in 
the top of the generator, it is forced by the low air 
pressure through a nozzle directed over the blade. The 
rotor is a one-bladed device driven by a variable speed 
electric motor. The speed of the blade and tunnel as well 
as control of the visualization devices is centrally located 
near the observation window. A strobe light with a variable 
flasher allowed the blade to be viewed once every revolution. 
Since the vortices are being shed at essentially the same 
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frequency as the blade rotates, they also appear to remain 
stationary in the free stream. Pictures of the flow were 
taken with a Polaroid camera using ASA 10,000 film. 

Figure 2 shows a schematic of the tunnel. Figure 1 is a 
picture of the symmetrical cross-section blade mounted 
horizontally. The blade is counterbalanced to insure rotor 
vibrational stability. 




Figure 1. Rotor Used in Horizontal 

Segments of Visualization. 
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Visualization Tunnel. 



III. THE GLAUERT EMPIRICAL CURVE 



In 1925 C.N.H. Lock and others proposed to augment the 
"vortex theory of airscrews" with empirical data. Since 
this theory breaks down for a range of working conditions 
between zero speed of advance and some negative speed of 
advance, a two-bladed airscrew was tested at angles of 
attack below the stall. The data from this airscrew was 
used to fill the gap in the range where the theory did not 
apply. In 1926 H. Glauert extended Lock's analysis to 
cover the complete range of working conditions. The data 
which eventually evolved is commonly referred to as Glauert's 
empirical curve. This curve is expr<=^‘=sed in ■f'ermc; of hwo 
non-dimensional coefficients, F and f, defined by the 
following relationships: 

= 4ttp (V^-v) (1) 

3— = 4ttp V f (2) 

dr ^ o 

is the speed of ascent or descent, v is the induced 
velocity at the blade element, and r is the distance of the 
element from the rotor axis. The results of these two 
investigations are given in Figure 3. 

The helicopter rotor in some instances acts as a pro- 
peller, notably in the hover and in vertical ascent or 
descent. The conventional anti-torque or tail rotor is 
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subjected to the same type of flow when the helicopter is 
in sideward flight. Therefore the empirical data is of 
interest in the study of helicopter aerodynamics. In the 
following pages the three working states of the rotor will 
be discussed and flow visualization pictures will be used 
to substantiate the discussion. 

A. THE NORMAL PROPELLER STATE 

The normal propeller state is the region on Glauert's 
empirical curve (Fig. 1) which starts at a point where 




Figure 3. Glauert's Empirical Curve. 
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1/f = 0 and where 1/F has one of several values. According 
to Glauert*s curve 1/F = 2.0; by the momentum theory 1/F = 
1.0. In practice 1/F is between these two extremes, nor- 
mally around the value of 4/3. The normal propeller state 
extends to the right as values of 1/F and 1/f increase. 

This state can be represented as shown in Figure 4 , T being 
the thrust of the rotor and u the flow through the disk. 
Initially, in the hover mode, = 0 and u = v. As power is 
increased, the helicopter ascends with a speed V^. The 
flow through the rotor is now + v and, from the momentum 
theory, is + 2v downstream. As the speed of ascent 

Vq+v 

increases , the value of u/V 1 since u/V = = 1+v/V . 

' o ' o Vq ' o 

In the limiting case of u/V^ = 1, 1/f = 1/F = Therefore 

■f — T? =: — rS — T7 

^ ^ v^. 




Figure 4. Diagram of Normal Propeller State. 

The photographs in Figure 3 were taken with the rotor 
arranged to simulate the normal propeller state. The rotor 
axis is mounted parallel to the free stream which flows from 
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left to right. The blade is rotating clockwise when viewed 
along the free stream vector. The thrust and induced 
velocity vectors are in the same direction as depicted in 
Figure 4. The tip vortices are obvious as is the vena con- 
tracta. The decrease in slipstream area downstream indicates 
an increase in velocity. In Figure 5-b the start of cyclonic 
tails which originate in the eye of the vortex are plain. 
Although not obvious in the second picture, these tails or 
trailing vortex lines extend downstream of the rotor in a 
sprial pattern, the visible series of tip vortices being 
simply the result of vortex flow around each spiral vortex 
line . 




a. b. 

Figure 5. Normal Propeller State. 
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1. The Hover 



As discussed previously, one of the limiting cases 
of the normal propeller state is the point where 1/f = 0 
and 1/F = 2. This is the hover mode and presents a severe 
challenge to the aerodynamicist in his attempt to predict 
rotor performance. The induced velocity is extremely 
important in this region because in a hover the induced 
power is approximately two thirds of the hover power for a 
helicopter with a conventional figure of merit. Previous 
attempts to calculate the rotor performance have not been 
very accurate because the vorticity from a blade is not 
swept away from the disk and strongly influences flow on 
the preceding blade. As the number of blades increases. 





Figure 6 . 



The Hover. 



b. 
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this effect becomes more pronounced. Figure 6 shows the 
blade mounted horizontally, with = 0 to simulate the 
hover mode. In the pictures the trail of vorticity is 
quite evident and is still in close proximity to the plane 
of the disk, even though the blade is past the beam posi- 
tion. With a four-bladed rotor, the next blade would be 
strongly affected by this vorticity. The pictures also 
show the vortices extending to the periphery of the disk 
and indicates the vena contracta under the rotor. 

B. THE VORTEX RING STATE 

If the helicopter is in a hover and power is decreased, 
it begins to descend under power. This is the vortex ring 
state and is represented empirically on Glauert's curve as 
the branch extending from the hover point, which has pre- 
viously been discussed, to the so-called point of ideal 
autorotation, that is, 1/f == 2.0 and 1/F = 0. Schematically 
it is represented by the following diagram: 



V 




V 



Figure 7. Diagram of Vortex Ring State. 
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In this case the thrust vector is in the same direction 
as the relative velocity in a descent. Initially, as power 
is reduced, the induced velocity v is much larger than the 
sinking speed and is negative and large. As the 

sinking speed increases, u/V^ becomes smaller and the 
induced velocity does not travel as far upstream before it 
is dissipated by the free stream. With further increase in 
V^, the tip vortex flow is driven downstream (with respect 
to V^) and is then drawn into the rotor flow. This recir- 
culatory phenomenon starts at the rotor tips and moves 
inboard. At some sinking speed, V^, the induced velocity 
can not escape from the plane of the rotor disk and the 
whole disk becomes embraced in turbulent flow. 

in Figure 8 the nuglc of attac]v has been modified to 
allow the test section to simulate a rotor in the vortex 
ring state. The free stream moves from left to right and, 
this time, the thrust vector points in the same direction. 
Figure 8-a simulates the rotor at a very low sinking speed. 
The tip vortices are clearly seen. The cloudy background 
indicates that the induced velocity has been dissipated 
and the flow is downstream outboard of the rotor disk. 

Figure 8-b shows the critical part of the vortex ring state. 
The blade is enveloped in a turbulent flow, the pattern of 
which is fluctuating periodically. The phenomenon results 
in a loss of thrust and severe vibration of the rotor. 
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a. b. 

Figure 8. Vortex Ring State. 



C. THE WINDMILL BRAKE STATE 

The last of the working states is the windmill brake 

state and occurs when the blades are extracting energy from 

the air. It is the upper branch of Glauert's curve from 

the ideal autorotation point to 1/f = 1/F = «* . Now 

u = V - V and u/V = 1 - ^ . As V increases, u/V 1. 

o 

Schematically this state is shown below. 




Figure 9. Diagram of Windmill Brake State. 
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Figure 10 is an extention of the previous pair of 
photographs except that the tunnel speed is increased above 
the value which caused the blade to be in the critical part 
of the vortex ring state. Both pictures clearly show that 
tip vortices are again formed and that the area downstream 
is greater than that upstream or the converse of the stream 
tube embracing the disk in the normal propeller state. 




a. b. 

Figure 10. Windmill Brake State. 
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IV. MANGLER'S DISTRIBUTION 



Glauert's assumption that the induced velocity is 
distributed uniformly across the disk is fairly restrictive. 
In 1948 Mangier calculated an induced velocity distribution 
which was more general. He assumed that there was an infi- 
nite number of blades, allowing the rotor system to be 
replaced by a circular disk. He also assumed that the tip 
speed ratio was large. This rules out investigation of 
the region around the hover mode, which has already been 
pointed out as a critical area. Mangier selected a plaus- 
ible loading over the disk (Fig. 11) and calculated the 
correspondina distribution of u. 




Figure 11. Non-dimensional Disk Loading. 
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The results are shown in Reference 4 for different values 
of disk incidence. The most significant deviation from 
Glauert's original assumption is that, since the disk is 
producing lift, there is circulation around the disk. Thus 
there is an upwash in the front of the disk. This phenomenon 
is shown clearly in Figure 12, there being an upwash in the 
front of the disk followed by a downwash. Notice also the 
vortices in the downwash througli which the blades are 
passing. 
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b. 



Figure 12. Mangier 's Distribution. 
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V. CALCULATION OF THE LINEAR STEADY STATE VELOCITY 
IN THE VORTEX RING STATE 

With the aid of smoke flow visualization, it is now 
possible to chose a plausible mathematical model in an 
attempt to calculate the induced velocity in the vortex 
ring state. These velocities will be determined by con- 
sidering the data from Glauert's empirical curve and 
applying a method of analysis similar to that used by 
Bennett in Reference 1. 

Consider Glauert's definition for the thrust 
coefficients : 

= 4-:rrp’.i^F (1) 

§ = ( 2 ) 

u is the resultant flow through the elemental annulus and 
is expressed as V^ - v in an area of upwash and v - V^ in 
an area of downwash , both expressions being positive. To 
calculate the induced velocity, a new thrust coefficient 
F will be defined: 

^ = 4TTrpv^F (3) 

Dividing equation (1) by equation (2) yields 
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Dividing equation (1) by equation (4) gives 




Substituting equation (4) into equation (5 ) , 




( 5 ) 



( 6 ) 



Re-writing equation (4) and expanding u, noting that 

(Vo - = ('^ - v„)^ 



I = (1 - 2v/V^ + (7) 

or 

1=1-2 Vf/F + f/F (8) 

Solving for 1/F 




F 



For the range of values for 1/f and 1/F in the vortex ring 
state, the value of 1/F can be determined. 



1/f 


1/F 


1/F 


0.0 


2.0 


2.0 


.5 


1.0 


2.91 


.6 


.75 


2.69 


. 8 


.50 


2.56 


00 

o 

1 — 1 


.25 


2.37 


2.0 


0.0 


2.0 
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Now Glauert's curve becomes three dimensional. The projec- 
tions are shown in Figure 13. The relationship between 1/F 
and 1/f is approximated by a fourth degree least squares 
polynomial fit. 

1 1^ 1^ l2 1 

= a(^) + b(j) + c(p + d(p + e (10) 




Figure 13. "Three Dimensional Glauert's Curve." 
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where a = -2.96 



b = 12.96 
c = -18.76 
d = 9.68 
e = 1.38 



Solving for f , recalling that ~ 

V F 
o 



f = i 

e 









, V ^ V ^ 

b O a. C O . 

^ - d 

r V V 



( 11 ) 



Substituting into equation (2) and multiplying the terms in 
2 

the brackets by V 



dT _ 47Tpr 
dr e 



■ 2 . b c 2 

+ ?3 — - i 2 -r + j — -<^^0 

r V F V F V 



( 12 ) 



Consider an expression for the thrust in terms of airfoil 
and dynamic pressure considerations. 



i= C,BCp(cor)2 1 



(13) 



Consider an elemental blade section, as in Figure 14. 

C^ = Aa = A(0+4>) (14) 

V^-v 

(15) 



Vq-v 

Co = A(6+ — — ) 
Z cor 



Substituting equation (15) into equation (13) and equating ' 
the result to equation (12) , recalling that 

2 

Raco rAe 



V -V V 


00 

> 

( 

L_ 


> 

> 






+ -§_ _2 O- 


o c o 


dV ^ 
o 


Kj ^ 1 — 

(Or j 


6 =2 


4 -F 2 " 


F V F 


V F V 


(16) 
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Figure 14. Elemental Blade Section. 
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c,wRxe + c,v - c, v = 

1 1 o 1 



2 a o b o . c o 2 

^ ^=3 — -=2 — •"= — 

F V F V F V 

(18) 

V ^ V ^ V 

2 ^ a o b o ^ c o 2 

V + — a ^ — r + T- ~ 

F^ F^ F 

(19) " 



From equation (4) 



m 

CN 

II 

P4 

CM 


(20) 


(V ^ - aV V + v^)F = V 
o o o 


(21) 


F-2^F + ^= f 


(22) 
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(23) 
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where 



, 2FF 

A = 33 

FF+fF-fF 

is the constant speed of ascent or descent, and v is 
the average value of induced velocity. Substituting 
equation (23) into equation (18) 

C^F^U3R0X+ C^F^vA- C^Fv = 

v^F^ - aA^v^ - bFv^A^ - CF^v^a"^ - dF^v^A^ (24) 

Finally 

KiV^ + K 2 V + = 0 " (25) 

where 

= (F^+ 2.96A^- 12.96FA®+ 18.76F^A^ - 9.68F^A^) 

K„ = Cj^F^(l-A) 

Kg = -CgF^uex 

Knowing the blade characteristics, the location on the 
three-dimensional thrust curve, and the non-dimensional 
distance along the blade, the induced velocity may be cal- 
culated. Some typical values are given below in tabular 
form and graphically in Figure 15. 



l/f 


1/F 


1/F 


A 


^average 


V 

0 


comments 


0 


2.0 


2.0 


0 


24.5 


0 


hover 


.5 


1.0 


2.91 


.42 


36.3 


15.1 


condition one Fig. 15 


.6 


.75 


2.69 


.47 


36.3 


17.2 




.8 


.50 


2.56 


.56 


33.1 


18.5 




1.08 


.25 


2.37 


.68 


30.8 


20.8 


condition ta*7o Fig. 15 


2.0 


0 


2.0 


1.0 


29.7 


29.7 


"ideal autorotation" 
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Figure 15. Flow Through Rotor for Various Points in the 
Vortex Ring State. 
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VI. CONCLUSIONS 



Flow visualization is an extremely useful tool which 
should be further exploited now that the facilities are 
available. More sophisticated methods of photography could 
be employed to record the vortices. Velocity measurements 
could be made in the wake with a device such as a hot-wire 
anemometer. These measurements, coupled with the photo- 
graphs, could help clarify some questions about the complex 
flow in the rotor wake, as well as other flow problems. 

Tip vortices exist in all the working states of the 
rotor just as an airplane wing has trailing vortices in all 
regimes of flight. i’ll tlic= hover mode tlio vortices Tint- 

escape from the plane of the rotor before the next blade 
comes along. The flow over the blade is, therefore, 
extremely complicated and rotor performance is difficult 
to calculate. 

At some point in the vortex ring state, the steady state 
velocity of descent prevents the induced vortices from 
leaving the blade. The whole rotor is enveloped in vortex 
rings which are symmetric about the axis. This condition 
is unstable and these vortices are shed in a non-uniform 
manner. 

Mangier 's distribution is shown to be at least quali- 
tatively correct. There is an upwash in front of the disk 
and the blade cuts through its own vorticity closer to 
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the hub. This blade-vortex interaction could cause unde- 



sirable vibrational effects. 

The induced velocity distribution for several points in 

the vortex ring state was calculated. One of these points, 

"ideal autorotation," is not an empirical point, as is the 

others, but a theoretical point from Reference 3. The 

values obtained appear to be optimistic, that is, half the 

blade in upwash and half in downwash. As reported in 

Reference 1 the value where V = v is much closer to the 

o 

tip. This calculation takes into consideration the profile 
drag of the blade. As disk loadings become higher, the 
figure of merit approaches one, the profile drag will 
become less important and the distribution will approach 
rhe idecil ca&c: reported in this paper. 
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